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ELECTRO-ARCO 


LIMITADA 
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POSTOS E ACESSORIOS 
PARA A SOLDADURA ELECTRICA 


MATERIAL APROVADO PELO 
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RUA DO BOLHÃO, 216- TEL 24277 
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Exija expressamente Papel Sihl 
Papeteries Zuricoises sur Sihl Zurich e 
Ma im 


Papel de desenho 
mate e acetinado 
granitado + 

para aguarelas, 


Papel de desenho 
acetinado. 
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ASSISTENT 
Sidi Papel de desenho 
DD a it mate, 

3 Papel vegetal 
Utoplex mate e acetinado, 
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COMPANHIA CIMENTO TEJO 
(Alhandra) 
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Assente sobre 407 estacas Franki 


UNDACOES FRANKI 
pittlcide onstrutora ortuguesa, ds, 


Praça do Areeiro, n.º 9, 4.º Esq. — Lisboa — Tel. 726061 e 729911 
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MOBILIÁRIO 
ARTALICO 


EM TODOS OS GÉNEROS 


Instalações completas para : 


Mobiliário moderno CLINICAS 
ROS PI-TAIS 


para 
SANATÓORIOS 
ESCRITÓRIOS CINEMAS 
LS COL AS EO UE SS 
BIBLIOTECAS ESPLANADAS 


SALAS DE EXPOSIÇÕES : 


Rua Febo Moniz, 2a 20 — Praça dos Restauradores, 49 a 517 


TEÊCNICA — II 


PIMENTEL & 
CASCAMRHO 
Quase 40 Anos 


AO SERVIÇO DA 
ENGENHARIA PORTUGUESA 
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(GRANDE SORTIDO 

DE INSTRUMENTOS DE 

PRECISÃO, CIENTÍFICOS 
E INDUSTRIAIS 


RUA vas PORTAS ve SANTO ANTÃO, 75 
LISBOA 
Te.er.: 24314-242 48 + TeLec.: TECNA 


SOCIEDADE | 
PORTUGUESA 
sk 


POSTES DE BETÃO PARA 
LINHAS ELÉCTRICAS 
E AZULEJOS 


MOSAICOS 


EMPRESA ELECTRO CERÂMICA—S.A. R.L. 
Candal — Gaia 


ISOLADORES DE CADEIA 


Para suspensão 
Tipo ACS 25 


COLUNAS PARA 
ILUMINAÇÃO PÚBLICA 


MANILHAS PARA ESGOTOS 


2 


RUA D. ESTEFÂNIA, 94-A 
TELEF. 47812-50129 


LISBOA 


Este isolador é igual ao tipo I. B. S. 99 da firma inglesa Taylor 
Tunnicliff & Co, Ltd., por acordo com a qual o fabricamos 


Peso aproximado: 
4,3 kg cada elemento 


Materiais 


Porcelana vidrada em castanho. 

Campânulas de ferro fundido maleável, galvanizadas por imersão a quente. 
Hastes de aço macio galvanizadas por imersão a quente. 

Molas de fixação de bronze fosforoso. 


Sede — Largo Barão de Quintela, 3 — LISBOA 
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Para um 
serviço permanente sem perturbações 
são ideais as 


pritadeiras de martelos DE ROLL, 
com duplo rotor, 


de construção moderna 
e sólida 
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Contadores de água fria para casas de habitação, diâmetro de passagem de |3 a 40 mm | 


EXPORTAÇÃO: 


ME TRIMFHENX 
4 Sociedade Húngara para o Comércio Externo de Instrumentos 


BUDAPE 31 


SREM p 


Correspondência: Budapest 62 Caixa Postal 202 — Húngria 


| Telegramas: Instrument Budapest! 


TÉCNICA— V 


MERCEDES-BENZ 


MOTORES INDUSTRIAIS 
DIESEL 
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MOTORES PARA TODAS AS APLICAÇÕES INDUSTRIAIS 

FAMOSOS EM TODO O MUNDO PELA SUA GRANDE ROBUSTÊS. ECONOMIA 

DE CONSUMO E FACILIDADE DE MANUTENÇÃO 
MODELOS DESDE 15 A 2.500 H., P. 
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MODELO MB 841 B ........ 25 H. P.— 1.000 R. P. M. 
MODELO M202 B ........ SS H. P.— 1.200 R. P.M. 3 
MODELO M 203 B ........ 90 H. P.— 1.200 R. P.M. = 
MODELO M 204 B ........ 120 H.P.- 1200R.P.M. 


OUTROS MODELOS PARA ENTREGA MUITO RÁPIDA 


MOTORES A 4 TEMPOS - ARREFECIMENTO POR AGUA 
ARRANQUE A AR COMPRIMIDO OU ELÉCTRICO 


REPRESENTANTES CG. SANTOS LDA. Divisão MARÍTIMA E TÉCNICA 


29, AV. DA LIBERDADE, 41 - LISBOA 
160, R. DE S.TA CATARINA. 168-PORTO 


TÉCNICA — VI 


ISOLAMENTO 
TÉRMICO DAS 
CONSTRUÇÕES 


PAVIMENTOS 


e 


TECTOS 


* 


PAREDES 
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CÂMARAS 
FRIGORÍFICAS 
a 
ISOLAMENTOS 
INDUSTRIAIS 


Terraços e pavimentos isolados com CELULITE no 
edifício da Comp.' de Seguros Fidelidade — P. das Águas Livres — Lisboa 


SETH, LDA. + RUA FILIPE FOLQUE, 10-1.º + TEL. 7301 56 — LISBOA 


Para impermeabilizar terraços, 
paredes, fundações, etc. 


COLAGEM DE 
TACOS DE MADEIRA 


TELEF. 730156  TELEG. EPALDA — LISBOA 


PL RTICA EM SAMOALAS EMA 
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MOTORES DIESEL 
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EQUIPAMENTO ELÉCTRICO 


Uma organização ao serviço 


das empresas armadoras 


ENGLISH ELÉCTRICA DE PORTUGAL, LIMITADA 
Rua Cascais, 47 (Alcântara) — LISBOA 3 
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O EMPREGO DO CIMENTO BRANCO... | 


permite acabamentos 

mais perfeitos, mais 

duradouros e muito 
mais económicos 


Estude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


Consulte os distribuidores gerais 


No Sul: SCIAL — T. do Corpo Santo, 15 — Telef. 2 0464 — Lisboa 
No Norte: SCIAL — R. do Bonjardim, 205 — Telef. 25779 — Porto 


mete lata tato Pintura ideal para paredes de prédios e outras estruturas. 
Po a O" 

elessiececa se Cores atraentes, fixas e duradouras. 
Cotefotadado? Usa-se sôbre rebocos, tijolos, betão, 
ou 00. 


pie la lalala Tel fibrocemento, pedra natural, etc. 
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DISTRIBUIDORES EM PORTUGAL pe tasera vera tenaracatatata tos 


HENRIQUES prspriiites 
& CASTRO, Za 

AVENIDA CONDE DE VALBOM. 96 LISBOA 
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MISTURA FRIA = Produz acima de 50 ton. | hora 
com o misturador 840-B. 


MISTURA AGREGADOS MÚLTIPLOS, DE ALTA QUALIDADE — Com rendimento 


MISTURAS QUENTES INTERMÉDIAS — Produz 
até 45 ton, / hora, utilizando o misturadar e secador, 


Niva Intalação do mistura dh 
befuminosos, marca BARBER 
GREENE, modelo 840-B 


Escolha a instalação de Mistura que mais 
o satisfará. — Por baixo custo inicial, po- 
de adquirir a nova instalação de mistura 
de betuminosos, marca Barber Greene, 
modelo 840-B que apresenta as caracte- 
rísticas próprias das grandes instalações. 
O misturador isolado produz bases esta- 
bilizadas e misturas betuminosas rias. 
Adicionando um secador, um seleccio- 
nador e um retentor de poeiras (ciclone) 
ao misturador, poderão produzir-se mis- 
turas betuminosas de auaisquers caracte- 
rísticas, 

O novo e mais comprido misturador, de 
veio duplo, propurciona uma mistura mais 
perfeita e uma maior capacidade. 


À tremonha de descarga do misturador 
é accionada hidrâulicamente, tacto que 
evita a segregação da mistura, À posição 
alta da tremonha de descarga permite que 
as camionetas parem para a carga, mesmo 
por baixo da tremonha. 

Facilmente móveis e montáveis, — Todas 
as unidades são transportadas sobre semi- 
«reboques equipados com rodas de pneus, 
facilmente rebocáveis por camionetas. 
Todas as unidades dispõem de suportes 
ajustáveis, não necessitando, pois, de 
bases, 


Para informações detalhadas dirija-se ao representante exclusivo. 


Borber-Greene Americos, Int, Avrora, Illinois, U.S.A, 
Borber-Greene Overseas, Inc, Auroro, Illinois, U.S.A. 


REPRESENTANTE EXCLUSIVO: 
SOCIEDADE DE MECANIZAÇÃO 
INDUSTRIAL E AGRÍCOLA, S.A. R.L. 


Avenida Padre Manuel da Nóbrega, 8 
SBOA 


SMEIA| 
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Barber-Greene é 


Barber-Greene Canodo, Ltd, Toronto, Ontario, Conoda 
Barber-Greene Olding & Co., Ltd, Hatfisld. Herts, Englond 


Ad No. 55-2034-B 
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ESCAVADORAS 


RUSTON-BUCYRUS, LTD. 
E 


BUCTIRUS ERIEGCO: 


Modelos de 3/8 a 36 jardas cúbicas 


TER CITE 
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As escavadoras de reputação mundial mais utilizadas no território português 


Stock permanente de peças sobressalentes 
Oficina de reparações 

Técnicos e mecânicos especializados 
Ensino e treino de operadores 


Representantes exclusivos para Portugal e Ultramar 


MONTEIRO GOMES, LIMITADA 


Rua Cascais, 47 (Alcântara) 


LISBOA 
LUANDA LOURENÇO MARQUES 
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Li» grande marca italiana 


Representantes em Portugal 


ELELCTR 


EMPRESA TECNICA DE EQUIPAMENTOS ELECTRICOS, S.A.R.L. 


Rua Sá da Bandeira, 766-1.' E. 


Rua Rodrigo da Fonseca, 103, rc 


Tel. 24818 - PORTO 


LISBOA — Tel. 686072 
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PINÇAS VOLT — AMPERIMÉTRICAS 


7 ESCALAS 


EM AMPÉRES EM VOLTS 


O-10 AMPS. O— 150 V 
O—- 25 , O- 600 » 
O — 100 
O- 250 
O — 1000 


PINÇA 
WATTIMÉTRICA 
7 ESCALAS 
O--3 KW 

0o—6 
0O—12 
0o—30 
Oo -—- 60 
O — 120 
O — 300» 


REPRESENTANTES: 


DIVISÃO MARÍTIMA E TÉCNICA 
C. SANTOS LDA. T. DA GLÓRIA, 17-LISBOA 


29, AVENIDA DA LIBERDADE 41 — LISBOA 160, R. STA. CATARINA, 168-—-PORTO 
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COMPANHIA PORTUGUESA 
DE FORNOS ELÉCTRICOS 


PRODUTORA DE 


FERREGUSA 
FERRO LCIGAS 
FERRO-MANGANÉÊS 


SÍLICO- MANGANÊS 
FERRO-SILÍCIOS 


15% - 25% - 45% - 75º - 90h 


CARBONETO DE CALCIO 
CIANAMIDA CÁLCICA 
PASTA PARA ELECTRODOS 


E EM MONTAGEM 


FÁBRICAS SEDE 
CANAS DE SENHORIM Largo de S. Carlos, 4 -2.º — LISBOA 
Telefone 67222 Telefones 25343-29608-368989 


AÇO TOR 


O AÇO DE ALTA RESISTÊNCIA 
PARA BETÃO ARMADO COM 
MAIOR VENDA NO MUNDO 


TENSÃO DE SEGURANÇA 2,400 Kg /cm 2 
SEM GANCHOS 


Especifique 
AÇO TOR 
com o símbolo 


) 


SOABAL - SOCIEDADE DE AÇOS PARA BETÃO ARMADO, LDA. 
RUA JOAQUIM BONIFÁCIO, 2, 1.º — TELEF. 405 66 — LISBOA 
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Director: JORGE GONÇALVES CALADO 
ADMINISTRADOR: MIGUEL FALCÃO 
secretário: ANTÔNIO ALVES MACEDO 


Ano XXXVI-N: 314 Julho de 1961 


C. D. U. 548 [378.962] 


CALCULO NUMERICO: Metodologia 


PELO ENG.o MANUEL FRIAS DE ALMEIDA E SÁ 


C. E.M. A.E.E.N., M.S.1,A.M. 


O Cálculo Numérico, oriundo das técnicas dos pequenos calculadores mecânicos e de 
adaptações ou criações esporádicas de certos modelos matemáticos, só nos últimos anos vem 
encontrando a sua estruturação específica neste crescendo incessante de quantificação do mundo 
material. Esta quantificação, porém, é dependente do conhecimento da estrutura física do 
mundo material, ainda incipiente em tantos domínios fundamentais da Tecnologia por falta 
de teorização do estado físico real e/ou por insuficiência de experimentação planificada (muitas 
vezes devida à dificuldade intrínseca da definição de certos processos materiais reais e corre- 
lativo deficiente controle dos mesmos). Esta problemática do quantitativo informa a Metodo- 
logia dos cursos, particularmente no Instituto. Foi, neste aspecto, que surgiram, à margem 
das lições da 51.4 Cadeira, os apontamentos e observações que seguem: 


1. O Cálculo Numérico é, por excelência, a tecnologia do quantitativo, devendo ser compreendido 
à luz da evolução científica e tecnológica das últimas décadas, em que se passou, um tanto brusca- 
mente, dum mundo em que a compreensão era um fim em si para uma atmosfera em que a compreensão 
é considerada, principalmente, promotora de acção de influência sobre o meio e/ou de controle deste. 

A Matemática, tendo evoluído por solicitações ou caminhos diversos, encontrou, na primeira 
grande ciência de quantiticação rigorosa, experientalmente confirmada — a Mecânica Celeste — motiva- 
ção para um desenvolvimento extraordinário da Análise e objectivo para um primeiro surto sistemá- 
tico de técnicas de cálculo. A Teoria da Gravitação, de Newton deve muito ou deve tudo ao estudo 
numérico das órbitas dos planetas, a que Kepler dedicou a sua vida de estudo, mas poucos apreciarão 
hoje quanto a criação dos logaritmos neperianos ajudou os trabalhos do astrónomo. É significativa 
a dedicatória que faz Kepler nas suas Ephemeris «ao ilustre barão Napier», de quem o pensador 
David Hume diria, mais tarde, que «the title of Great Man is more justly due than to any other 
whom Scotland ever produced». 

Só no fim do século passado se deu, porém, uma autêntica sistematização da Matemática. 
Lembremos, apenas, a unificação e sistematização da Análise por Cauchy e o célebre programa de 
Erlangen, de Felix Klein, que mostrou a íntima associação das Geometrias com a Álgebra. Própria- 
mente, só já no segundo quartel deste século — e depois de um trabalho exaustivo de classificação, 
unificação, arrumação e exclusão — encontrou a Matemática o seu objecto específico — o quantitativo, 
axiomatizado, não contraditório — e o seu método próprio — o método algébrico. Ficou, assim, definiti- 
vamente estruturada a Matemática como linguagem do quantitativo, como ciência independente, essencial- 
mente não experimental, com objecto e método próprios. Já J. W. Gibbs, numa visão genial, deixara 
entrevisto, nos seus escritos precursores da moderna interpretação estocástica da fenomenologia, 
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este carácter essencial da Matemática como linguagem do quantitativo. Tecnicamente, esta evoluiu 
segundo dois ramos ou especializações profissionais — Álgebra e Topologia. 


2. O Cálculo Numérico é Matemática apenas no aspecto formal. Como tecnologia do quanti- 
tativo tem que respeitar e empregar constantemente a linguagem do quantitativo, que é a Matemática. 
A Física, análogamente, enquanto quantitativa tem que utilizar a linguagem matemática. As equa- 
ções de Maxwell, do electromagnetismo, têm uma expressão matemática porque traduzem interde- 
pendências quantificadas de variáveis físicas mas a sua essência é não matemática. A equação funda- 
mental da Mecânica Clássica tem uma expressão matemática mas é, essencialmente, uma síntese 
experimental QU), 

Num mundo em que só nas últimas décadas se deu uma quantificação apreciável da experi- 
mentação e das Tecnologias, era natural que se enquadrasse sob a designação geral de Matemática 
todos aqueles domínios de quantificação mais intensa, isto é, com linguagem quantitativa mais densa, 
Assim, ainda é hoje corrente apelidar de Matemática a Estatística, ciência que estuda o comporta- 
mento aleátório da fenomenologia material e que em nada contribuiu para a evolução daquela. Apenas 
se pode considerar matemática aplicada a teoria moderna da Probabilidade. Esta é, para a Estatística, 
simplesmente uma medida (medida da incerteza da generalização amostra — população) e, como tal, 
foi objecto de tratamento matemático aplicado, axiomático, como o integral definido. 

Ainda a respeito da confusão entre o formal e o fenomenológico é curioso citar a justifica- 
ção dada por P. G. Hodge para apresentar um tratamento formal da Plasticidade (2): «A palavra 
«plasticidade» significa muitas coisas para muita gente... e é considerada aqui um ramo das mate- 
máticas aplicadas». Ora o que é facto é que todos os que conhecem, por experimentação, a feno- 
menologia da Plasticidade sempre se entenderam sobre o sentido físico de tal termo. 


3. Metodologia é aproximação e introdução — aproximação porque foca os problemas sob uma pers- 
pectiva particular e introdução porque é apresentação a outrem de um novo domínio. A perspectiva 
metodológica do Cálculo Numérico evoluiu mais nas duas últimas décadas do que nos 3 séculos ante- 
riores, em que, desde os primeiros cálculos das efemérides, pelos astrônomos e analistas tiveram 
os cálculos muméricos uma evolução lenta e quase apenas subsidiária das matemáticas puras. 
As solicitações de um business contabilizado, ordenado mais por convenção do que por ciência, mas» 
aspirando a uma dinamização jamais conhecida e a uma produtividade vital, levaram, já nos fins do 
século passado, a um grande impulso da mecanização da aritmética contabilística, de que foi precur- 
sora a histórica máquina aritmética de Pascal. Criaram-se, assim, grandes empresas dedicadas ao 
fabrico de calculadoras mecânicas e ao estudo da mecanização do calculo aritmético. 

Há duas décadas, porém, o cúlculo automático e a automação industrial eram realidades incipientes, 
praticamente sem «peso» na estruturação da vida científica e na organização técnica. Hoje a progra- 
mação em termos de operações elementares — aritméticas ou unitárias — e o agrupamento delas em rotinas 
automatizáveis são o new look do Cálculo Numérico e da Tecnologia. 

Vencida a inércia mecânica dos primeiros grandes calculadores digitais pela substituição dos 
circuitos mecânicos por circuitos electrónicos, encontramo-nos, nas duas últimas décadas, em presença 
de uma corrida, pela parte das grandes empresas construtoras de calculadores, às maiores veloci- 
dades possíveis de operações elementares e à criação de programações para os cálculos comerciais 
e científicos o mais eficientes possíveis, dada a capacidade e a organização lógica das máquinas. 
Constituem-se, deste modo, as bibliotecas de programação, IBM, Bull, etc. 


4. A teoria da interpolação tem hoje um interesse crescente devido à necessidade de adaptar as 
estruturas continuas da Matemática, apoiadas sobre o conceito de número real a modelos discretos, 


(1) Não consideramos, é claro, certos pontos de vista críticos como os da corrente de Sir James Jeans da Física- 
“Matemática, nem podemos entrar aqui em certos aspectos da problemática actual da mesma, expostos pelo sr. prof. 
eng.” António Gião na sua lição do último Curso de Extensão Universitária, da Universidade Clássica de Lisboa. 

(2) in «Elasticity and Plasticity», Wiley, 1958. 
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construídos sobre o conceito de múmero inteiro e informados pela lógica da álgebra de Boole. 

Newton, que depois de Gregory, foi um dos criadores da teoria da interpolação, estabeleceu 
logo um dos esquemas mais gerais e completos de interpolação polinomial na sua teoria das diferenças 
divididas. Foi pena que Newton não tivesse sido tão feliz na escolha de notações como o foi no cálculo 
das fluxões. A sua notação não põe em relevo o carácter operacional do método das diferenças finitas. A fór- 
mula de Lagrange da interpolação, considerada um resultado fundamental em Álgebra, não é senão 
uma forma particular da formula de Newton. 

Um alto ponto da teoria da interpolação, no seu período clássico, foi atingido com a ideia de 
Gauss de escolher convenientemente as abcissas de interpolação de modo a obter-se a melhor apro- 
ximação para o integral. Foi uma ideia frutuosa, generalizada a intervalos infinitos e aplicada às 
funções mais importantes. Descobriram-se ou estabeleceram-se ligações desta teoria com a teoria das 
fracções contínuas e com a teoria dos polinómios ortogonais. A discussão dos valores próprios e das 
funções próprias levou às ligações recentes com a teoria das equações integrais e com estudos 
numéricos sobre a equação de Schroedinger. 

Um novo impulso na Teoria da interpolação foi dado por Tchebychev que introduziu a ideia 
da escolha dos coeficientes de interpolação de modo a obterem-se coeficientes iguais na fórmula 
integrada. Uma vez que estes coeficientes são os «pesos» com que os dados entram na fórmula, 
este ponto de vista estava de acordo com a, então muito em voga, teoria dos erros de observação. 

O período moderno da teoria da interpolação é caracterizado pela introdução dos pontos de 
vista da análise funcional, Radon, particularmente, estabeleceu uma fórmula muito geral para o 
termo residual («error term») da interpolação polinomial. Em muitos casos, este pode obter-se sob a 
forma monomial Ct, hm ff?) (L), sendo h o intervalo de interpolação e ff) (£) um valor inter- 
médio da derivada de ordem p da funcão f(x), em questão. Noutros casos é necessário recorrer a 
dois ou mais monómios. Ch. Blanc introduziu a ideia do estudo da «estrutura estatística» das fun- 
ções. Isto mostrou, entre outras coisas de interesse para o Cálculo Numérico, que a fórmula de 
integração de Simpson é, em geral, mais favorável do que as outras fórmulas de Newton-Cotes. 

Particular interesse têm hoje os métodos destinados a «acelerar» a convergência das iterações. 
É curioso notar que, no caso da inversão de um operador linear, obtém-se uma série 1 + K + 
+ Kº + Kº +...., de convergência linear; esta convergência pode passar a quadrática desde que 
se atenda ao produto de Euler (1+K) (1+Kº) (1+4+K*)... = RE 

Era corrente considerar-se como inconveniente para o cálculo toda a expressão analítica sob a 
forma de quociente cujo denominador fosse muito pequeno. Ora muitas derivadas aparecem sob 
esta forma. Este problema não existe, praticamente, para o cálculo digital moderno. Se se quiser 
obter um valor aproximado de tal quociente, por exemplo, com 10 decimais, basta tomar o nume- 
rador e o denominador com 20 decimais, isto é, trabalhar com precisão dupla. 

Metodológicamente é vantajoso—e começa a ser corrente — estudar a interpolação através de 
operadores, que, pondo em relevo a estrutura da interpolação, facilitam o domínio dos métodos e a 
adaptação destes ao cálculo com máquinas. Nos últimos anos, particularmente, sob a iniciativa de 
Znedek Kopal —- um grande nome da análise numérica —, têm-se desenvolvido métodos operatórios 
de integração e diferenciação numérica em termos de diferenças tabulares por um processo linear suscep- 
tível de notável precisão, sem recurso a diferenças finitas de ordem superior, As «operações» de diferen- 
ciação e integração são traduzidas em termos de operadores de deslocamento lateral e substituidas 
por aproximações em termos de fracções, racionais, com convergência rápida, A fim de se linearizarem 
as operações, as fracções são eliminadas por multiplicação. A aproximação racional do operador integral 
por diferenças centrais e diagonais permite uma integração sistemática das equações diferenciais ordinárias, 
Deste modo se consegue uma apreciável redução de tempo e trabalho de programação e cálculo. A maneira 
corrente de derivar as fórmulas de interpolação e integração numérica em temos de diferenças fini- 
tas é desenvolver as expressões E” ou log E, em que E é o operador diferença lateral, em potências 
crescentes dos operadores. Se, em particular, as expressões (1 + 4)” ou (1— v)”" são desenvolvidas 
por meio do teorema do binómio, obtêm-se as conhecidas fórmulas de Newton para a interpolação. 
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Se, analogamente, se se desenvolve log E em série de Taylor e se eleva à potência conveniente, 
todas as conhecidas fórmulas de integração e diferenciação se obtêm. Todas estas fórmulas possuem, 
contudo, um grande inconveniente: a convergência deixa muito a desejar e obrigam a empregar, por 
isso, várias diferenças de ordem superior. Em vez do desenvolvimento em série de Taylor, Znedek 
Kopal emprega uma aproximação por fracções cujos termos são polinómios nos operadores. A frac- 


ção continua do segundo membro de 
A 
log (1 + 4) = a 


LP SRA 
3 + “— Q2A 


4 


dá, por, sucessivos cortes, origem à sucessão. 
2Á 
| 
2 + À 
64 + 34º 
64644 4º 
R 60 À + 60 4º 4 11 4º 
E sacas ion 
60 + 904 + 364º + 34º 
A A? A3 As 
ARO NA RSS E Ta . Comparando esta sucessão com a série de Taylor, 
420 + 840 À 4- 540 4º 4-120 4º 4. 6 4+ 
de log (1 + 4) é fácil mostrar a superioridade dela sobre esta quanto à rapidez da convergência. Partindo 
a! 
da relação log (1 + 4) F=hf e fazendo log (144) =Rg = e . ; vem a relação aproximada 
"q; 
Ro(d)F-hf, que é não linear em À, Pods, contudo ser linearizada, suprimindo as fracções : 
Na (4) F-hD, (2) f. Vêm, as relações lineares 
2 dos=h (2f + 405) 
6ds + 3d =h (6fo +6 A, + 4) 


60 dy; + 60d, + 11d, =h(60fo + 90 4, + 36 4, + 3 4j,), etc. 


Ra = 


——— 


A tradução destas relações em termos do operador diferença central, ? (que mantém, é claro, 
a linearidade) permite obter fórmulas de integração numérica e, particularmente, um método de 
integração numérica das equações diferenciais ordinárias muito eficiente pelo relativo pequeno número de 
operações elementares em que se traduz a sua aplicação. 


9. A solução das equações integrais interessa não só às ciências puras mas, principalmente, à 
Tecnologia pois encontramos equações integrais na expressão de relações quantitativas em Plastici- 
dade (fenómenos com hereditariedade), em Hidroaerodinâmica (movimento dum fluído viscoso numa 
conduta, vibração duma asa pela tensão dum fluído incompressível), em Termodinâmica (distribuição 
de temperatura numa barra em regime não estacionário), em Óptica, em Oceanologia, etc. É inte- 
ressante saber como o cálculo numérico moderno veio permitir obter soluções expeditas para os problemas 
formulados em termos de equações integrais: 

Resolver as equações integrais é problemo matematicamente muito complicado, como é sabido 
do Cálculo Infinitesimal. Portanto, no Cálculo Numérico, é preciso limitar o tipo de equações a resol- 
ver pela máquina. A limitação diz respeito, particularmente, ao tipo de núcleo da equação integral, 


k (x,5): 


b oux 


2 (6) =F(x)— [k (xs) 2 (5) ds 


TEONICA 
582 


Como se sabe, b ou x, determinam os tipos Fredholm e Volterra, e v (x) = ou FO dá a espécie, 
1.2 ou 2.2 respectivamente. Toda a equação diferencial de ordem arbitrária e coeficientes constantes pode 
reduzir-se a uma equação integral de Volterra com um núcleo do tipo k (x— ss) e, por isso, as equa- 
ções com um núcleo deste tipo têm sido particularmente estudadas no cálculo numérico. Na deter- 
minação das características dinâmicas das soluções das equações integrais usam-se como dados iniciais 
dependências obtidas experimentalmente (1), Uma vez que as curvas experimentais são obtidas com 
erro de 5 a 10 “/y, este erro, vai reproduzir-se na solução numérica, que é, então, necessáriamente 
simplificada. . 

Os métodos usados para resolver as equações integrais são iterativos, do tipo simples e do 
tipo Gauss-Seidel. Uma equação de Fredholm de 2.2 espécie pode formular-se, então, assim, para O 
método de iteração simples 

(1) b (k—1) 


PO)=tf6)—1 | k(x 5)? (s)ds ; 
O 
e para o método de Gauss-Seidel : 
x=A 
(k) n (6) a (e=1) 
e t)=f6)—) [k(xs)? (s)ds—2 | k(xs)2 (s) ds, 
o Xi 


em que o índice superior entre parêntesis é o índice de iteração e À é uma quantidade arbitrária- 
mente pequena, em princípio. Como as fórmulas mostram, a diferença entre o método de Gauss- 
-Seidel e o método simples é que, naquele, ao calcular a função desconhecida, tomam-se os valores 
iniciais justamente calculados no passo iterativo em causa enquanto no método simples se empregam 
os valores de função incógnita obtidos no passo anterior. Por isso o método de Gauss-Seidel é mais 
conveniente. Particularmente, para equações de Fredholm de 1.2 espécie e com núcleo simétrico, de 
forma quadrática definida positiva, a convergência do método fica preestabelecida. 

Têm um grande interesse hoje as equações integrais de 1.2 espécie cujos núcleos são funções de 
autocorrelação e, portanto, funções simétricas, definidas positivas. Portanto : 


Db b 
[ [rx ol) o(s) de ds>o, 


onde 9 (x) é uma função arbitrária não idênticamente nula, Nestas condições, há, quando muito, umá 
solução para a equação Fredholm, pois, se houvesse duas, vw (x) e g2 (x), era 


b 
/ k (x, s) [9 (5) — ca (s)] ds =0 
b b dá 
e f f k (x, S) [24 (x) sas a (x) | [24 (5) e a (s)] ds da =: | 


o que era contra a condição anteriormente estabelecida. 

O facto de que uma equação de 1.2 espécie tenha uma só solução possível permite substituir 
o integral por uma soma o que conduz à solução dum sistema de equações algébricas cuja matriz é simé- 
trica e definida positiva. Suponhamos que a função w (x) tem a forma 


= tas 


wo (x) = 2 E 0 (x— xi) 


(1) V/ V. Solodovnikov, Vvedeniye v statisticheskuyu dinamiku sistem automaticheskogo upravleniya, Gos- 
tekhizdat, 1952, 
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em que 3(x) é uma função de Dirac e os É; são números arbitrários não simultâneamente nulos. 


Vem então: 
b b 


b b | 
f fr, 5) w (x) o(s) dx ds = 5) wo (x) dx ) k (x,5) S(s—si) ds = 


t a 


b 
n 


n 
E | o(x) dxk(x, si) = 2 2 E E | k (x, si) 9 (x— x;) dx = 
a ii A 


= [ím 


n nm 
2 Et k(x,s) 2h k (us) >0. 
j 


- [ás 
ma 


— bo 


No caso da função de autocorrelação, além de termos simétrica a matriz, k (xi, sj) = k (x;, si) 
é ainda k (x;,si) > 0. Então, o método de Gauss-Seidel é necessariamente convergente e podemos 
empregar a programação para os sistemas de equações sem qualquer condição particular. 


6. A introdução dos métodos da Estatística, particularmente nas comunicações e no controle 
automático, exige um cálculo intenso de uma característica fundamental dos processos estocásticos 
— a função de autocorrelação. Pode calcular-se em simuladores electrónicos e em calculadores digitais. 
No primeiro caso só é possível o cálculo para um número limitado de processos. O problema é radi- 
calmente resolvido pelo emprego de pequenos calculadores especializados no cálculo do valor médio 
M (x) e da autocorreloção R (7) dum processo estacionário x (t): 


+ oo T 
Mit [e lide = Hm a (di mem; 
T>% 2T). 


em que w(x) é a densidade de probabilidade do processo estacionário e T é o tempo de observação 
do processo x (t). 
O cálculo de função de autocorrelação é baseado no teorema ou princípio da ergodicidade: 


T + oo 
R (7) = Jim pe [x (t+-x(t)dt=M[x(t+-)x(D]= E (t+ T)xa(t)w(x1,x2;7) dxi dxs 
—T = Fe 


Na prática observamos de O a T, a informação sobre o processo é incompleta e, portanto, 
temos, aproximadamente : 


É 
E 
Ê t) dt 
m 7 )* 8 


T 
R (5) = Tjx6+9)x0 dt 


ê) 


É suficiente determinar R (x) entre O e T porque a função de autocorrelação é par. 

A função de autocorrelação R (7) é então calculada em calculadores contínuos com orgãos para 
produzirem uma desfasagermn 7 e para multiplicarem x (t) por x (t + 7) para vários valores discre- 
tos de 7 e, também, um dispositivo para integrarem o produto x (t + 7) x (t) durante o tempo T, de 
observação do processo. É preciso fixar o curso do processo num oscilograma ou num fonograma 
magnético para produzir a desfasagem. Trata-se, assim de calculadores intermédios entre os ana- 
lógicos e os digitais — os correladores. Entre estes há a salientar o fotoeléctrico, em que o processo esto- 
cástico é convertido em sinal eléctrico por meio de células fotoeléctricas. 
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7. Um dos aspectos mais representativos da quantificação do mundo material é, presente- 
mente, o da introdução dos métodos da Estatística em quase todos os domínios da experimentação 
e da previsão. Isto corresponde, por um lado, à necessidade de tomar decisões em face de um número 
limitado de dados numéricos e, por outro, a uma evolução na interpretação da fenomenologia 
física. A estruturação da Física —e, em particular, da Química — é essencialmente dedutiva e idea- 
lizada. Para o experimentador, para o tecnólogo, a Física é fundamentalmente a ciência das tendências 
das variáveis (físicas), sendo a Química uma Física com transformação de matéria. O experimentador 
e o tecnólogo, em especial, têm que colher dados, observar, controlar e prever não em condições 
ideais mas em situações reais. Têm que decidir sobre um complexo interdependente de multiplas 
variáveis físicas, com um número limitado de observações ou dados. 

Para controlar um processo, é essencial medir os vários valores que o caracterizam. Definamos 
a quantidade de informação relativa à leitura de uma variável, por exemplo, temperatura, pressão, tensão 
etc. Suponhamos que a quantidade varia entre os limites xa = Xmáx. — Xmin.. Mesmo para uma varia- 
ção contínua de x, o instrumento de medida fica insensível a uma variação inferior a + Ax, de 
acordo com a sua capacidade de resolução. Portanto só é possível obter um número discreto de va- 


Xd A = a % e a) 
lores de 4x: No Por consequência, a quantidade de informação associada a uma 
p X 


leitura é 


[= loga N = logs (+ 1) 
A x 


Como resultado dum processo industrial, deve obter-se um produto final satisfazendo a certas 
características técnicas. As características físicas do produto final, por cujas variações se pode julgar 
até que ponto foram satisfeitas condições técnicas, constituem a informação para o departamento do 
controle técnico. A selecção da comunicação essencial a partir de todo o agregado de características 
do produto é conduzida através da tecnologia específica da produção, que estabelece uma sequência 
no processamento da matéria-prima ou componentes primários. Os meios técnicos empregados para 
mudar as características da matéria-prima são analogicamente correspondentes ao transmissor dum 
sistema de comunicação. A matéria-prima ou os componentes primários, cujas propriedades físicas 
contêm a informação essencial que constituem o canal de comunicação. A parte do receptor no pro- 
cesso industrial é desempenhada pelos instrumentos que controlam a qualidade, a partir dos quais 
se sabe o produto final satisfaz os requisitos técnicos. Estes instrumentos consistem, basicamente, 
em duas partes — uma parte de medida e outra de filtro. O produto satisfatório corresponde a um 
sinal útil e o insatisfatório a uma interferência. Deste modo os instrumentos de controle são filtros que 
deixam passar as comunicações essenciais e eliminam os ruídos. 

Num sistema de comunicação nem o emissor nem o receptor tomam parte activa no funcionamento 
do sistema. Em particular o problema de saber como o emissor obtém a informação não se põe. 
O receptor tem um papel passivo e não interfere no processo de transmissão da comunicação mesmo 
que este seja insatisfatório. Por outro lado, num processo industrial controlado automaticamente 
tanto o emissor como o receptor desempenham um papel activo no processo da comunicação, que se 
traduz em acções em ciclo de retorno (feedback) entre o canal de comunicação e o emissor e também entre 
o recipiente da comunicação (isto é, o departamento de controle técnico), o transmissor e o receptor 
(isto é, o processo de fabricação). 

Há três teorias, à primeira vista muito distintas, para o estudo dos sinais acompanhados de 
interferências — sinais estocústicos —: a teoria da informação, de Shannon, a teoria de autocorrelação e filtração 
de Wiener e a teoria do teste de hipóteses. Estudos recentes evidenciam a sua interligação. A teoria de 
Wiener tem obtido um grande êxito nas comunicações mas tem dois defeitos essenciais: não conduz 
a uma optimização da filtração de acordo com os recursos actuais do Cálculo Numérico e parte de 
um conhecimento apriorístico, estatístico, do sinal e do ruído. Por isso, a teoria dos jogos, de von 
Neuman e Morgenstern e a teoria das decisões estatísticas, de Wald têm sido objecto de aplicação na 
«filtração» dos sinais. Na teoria dos jogos há dois adversários, um dos quais corresponde às condi- 
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ções exteriores a que o sistema com controle automático deve reagir; o outro adversário é o próprio 
sistema. Na teoria das decisões admite-se que o emissor pode gerar n sinais diferentes, um dos quais 
entra no sistema, acompanhado de ruído. O problema é isolar um dos sinais do complexo sinal- 
-ruído, o que é conseguido, teóricamente, através duma função de decisão. 

A teoria dos jogos está intimamente associada à programação linear que numéricamente se 
produz na optimização dum sistema de inequações com variáveis não negativas, para o que os cal- 
culadores numéricos digitais têm um programa apropriado. 


8. A experiência numérica é um novo tipo de experimentação que interessa à optimização em 
Tecnologia, particularmente naqueles domínios em que a estrutura quantitativa dos processos indus- 
triais se encontra bem estabelecida mas em que há necessidade de contar com a maior ou menor 
incerteza de certos coeficientes empíricos. A rapidez dos calculadores digitais veio permitir que 
em muitos casos a experiência em instalação piloto pudesse ser substituída ou, pelo menos, tornada 
mais eficiente com uma optimização prévia de condições. É disto exemplo flagrante o projecto de 
uma coluna de destilação. 

Os métodos do Cálculo Numérico são, presentemente, quase o único meio de estudo dos re- 
gimes de plasticidade em meios contínuos, por exemplo, em placas. Só em casos muito excepcionais 
elaborou a Análise Matemática modelos quantitativos adequados. Em estruturas discretas têm sido 
utilizados com êxito os métodos da programação linear e da programação quadrática. 

Os métodos do cálculo numérico têm tido um interesse particular no estudo do movimento 
dos satélites artificiais, cujo curto período de revolução (o Vanguard I, lançado em Março de 1958 
já deu mais voltas à Terra do que esta deu à volta do sol desde o início das civilizações!) inibe 
praticamente que se consiga estudar o seu movimento por integração das equações da Mecânica 
Celeste. O problema tem sido resolvido considerando a órbita do satélite como uma elipse instan- 
tânea, cujos elementos variam continuamente, de acordo com uma «função de perturbação», que é 
estudada numeéricamente. Já o mesmo tinha sucedido com a teoria do movimento da Lua, cujo movi- 
mento Newton, Euler, Laplace, Delaunay, Brown e muitos outros matemáticos e astrônomos não 
conseguiram explicar com toda a sua análise matemática. A teoria de Delaunay, por exemplo, incluia 
desenvolvimentos analíticos com centenas de termos periódicos. Por isso o estudo do movimento 
da Lua é hoje feito por métodos puramente numéricos que dão soluções exactas para um período 
de tempo limitado, embora não dêem conta da natureza dos factores intervenientes no movimento 
do planeta. 

Um dos tipos de experiência numérica mais interessantes e que tem encontrado muitas 
aplicações é o da simulação estatística, conhecida por método de Monte Carlo, em que se faz uma 
experiência numérica, «observando» se pontos dum espaço uni- ou multidimensional pertencem ou 
não a um domínio física ou matematicamente bem definido. Assim se têm resolvido complexos 
problemas de Investigação Operacional, de integração definida, de integração de equações às deri- 
vadas parciais e problemas de Química como o da determinação das funções de distribuição estatis- 
ticas de um conjunto de moléculas, por uma generalização da equação de Boltzmann. Um dos domí- 
nios de maior aplicação da experiência numérica pelos métodos de Monte Carlo é o da solução das 
equações de transporte nos reactores nucleares, para cujos problemas numéricos foram justamente cria- 
dos tais métodos. 


9. A organização dum programa actualizado de Cálculo Numérico, Mecânico e Gráfico para 
quem começa e em meio ano lectivo é um autêntico problema de escolha. O mais simples é subor- 
dinar o programa à Matemática, com todo um cortejo de teoremas e demonstrações. Mas, com isto, 
perder-se-á o essencial do Cálculo Numérico como tecnologia do quantitativo «You, mathematicians, 
know how to solve this problem but you can't actually do it» (1), Tendo presente a evolução da 
problemática do quantitativo, pareceu conveniente a elaboração do seguinte programa: 


(1) E, Milne, Numerical Calculus, 1949, pág. V, 
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MAGUE 


CONSTRUÇÕES METALO-MECÂNICAS MAGUE, L.?* 
(S. E. MONIZ DA MAIA & VAZ GUEDES, L.*) 


Alverca — Portugal 


Guindastes eléctricos de 5t/16m em funcionamento no Porto de Leixões 


PROJECT Q ET ADRICO DE 


GUINDASTES 
GRUAS. DERRICKS E GUINCHOS 
APARELHOS DE ELEVAÇÃO ESPECIAIS 
CONSTRUÇÕES METÁLICAS PESAD 


COMPORTAS | 
TURBINAS 


CONDUTAS 
FORÇADAS 


E BLINDAGENS 


DIVISÃO 


METALO-MECÂNICA EE 


COMPANHIA UNIÃO FABRIL 


TÉCNICA — XVHI 


AVENIDA 24 DE JULHO, 


170 —-LISBOA 


comporta 

da barragem 
do Maranhão 
construída pela 
companhia 
união fabril 


I — Elementos de Cálculo Mecânico (a fim de se conhecerem as potencialidades e programação 
dos calculadores) ; 
[I — Diferenças Finitas. Interpolação ; 
II — Equações às Diferenças e Equações Diferenciais; 
IV — Métodos de solução dos sistemas de equações lineares ; 
V— Cálculo numérico de séries e integrais ; 
VI — Complementos ; 

VII — Solução numérica de problemas de Engenharia (sistema com n graus de liberdade e ligação 
elásticas, resposta dum sistema não linear, distribuição de temperatura em placas con- 
dutoras, etc.); 

VIII — Elementos de Cálculo Gráfico. 


Procurou assim dar-se concretização à concepção do Cálculo Numérico como ciência autónoma 
tendo como objecto próprio o estudo dos sistemas discretos e como método fundamental o método das 
diferenças finitas. No estudo preliminar de Cálculo Mecânico encontra-se matéria de particular interesse 
pedagógico como os conceitos de circuito com retorno, programação e automação. E a muito mais 
poderá conduzir o estudo da organização lógica dos calculadores. O Cálculo Gráfico fica limitado a 
um estudo de Nomografia pois mais o não permite o tempo nem a unidade em curso. 


10. O mundo do Cálculo Numérico é todo o mundo da quantificação da matéria e da energia. 
Esta quantificação é essencialmente dependente do conhecimento da estrutura física do mundo mate- 
rial. Sem este conhecimento ou subordinação poderá cair-se no caminho perigoso da conjectura 
numérica, como faz por exemplo o sr. M. S. Klamkin (2) que, admitindo que um bife é homogêneo, 
isótropo e homotético chega — sem experimentação quantitativa — a uma lei quantitativa para o tempo 
de trabalho culinário e, não contente com a sua lei, generaliza-se para um bife a n dimensões! (que 
ele nunca deve ter comido); ou então como faz o sr. P. G. Hodge que, afirmando ninguém se enten- 
der sobre o estado plástico, resolve considerar este como um capítulo matemático (v. número 2). 

Como ÂAnteu que recebia novas forças ao tocar na terra, poderá a Matemática encontrar novas 
directrizes ao debruçar-se sobre os problemas numéricos actuais. O ensino do cálculo numérico deve 
ter em conta, porém, que o quantitativo se subordina ao sentido do qualitativo que é o da fenome- 
nologia física. E não pode esquecer que a metodologia tem as suas limitações próprias e exigências 
especificas, como aproximação e como introdução. 


(2) SIAM Review, Abril, 1961, pág. 167. 
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6.º Congresso Internacional de Preparação de Minérios 


Cannes, 26 de Maio a 2 de Junho de 1968 


Realizar-se-á em Cannes (Alpes Marítimos, França) de 26 de Maio a 2 de Junho de 1963, o 
6.º Congresso Internacional de Preparação de Minérios, que será acompanhado de visitas técnicas 
em França e noutros países mediterrânicos. 

O tema geral do Congresso será a Preparação de Minérios sob todos os seus aspectos, sendo dado 
especial projecção a 


— Concentração dos minérios de ferro, compreendendo instalação magnetisante, mas excluindo 
aglomeração e redução directa; 

— Medições, controle, regulação, automatização; 

— Hidrometalurgia. 


Podem ser apresentadas, como é óbvio, memórias sobre outros assuntos não mencionados nos 
três temas principais. 

Quando da escolha definitiva das comunicações, a comissão científica apoiar-se-á, antes de 
tudo, sobre o critério de valor científico e técnico. Mas, em caso de qualidade e interesse iguais, 
serão aceites de preferências, se for necessária uma escolha, as memórias tratando de um assunto 
referente aos temas principais. 

A comissão científica decidiu dar o maior valor a toda as considerações de ordem económica, 
que poderão ilustrar as comunicações, como sejam: custo de instalação, custo do tratamento, com- 
paração económica de diferentes processos concorrentes. 

Aos Investigadores e Engenheiros susceptíveis de apresentar qualquer comunicação a este 
Congresso, a comissão organizadora pede para estabelecerem, o mais cedo possível, contacto com o 
Secretariado do Congresso. 

Toda a correspondência deve ser endereçada a 


Sécrétariat Général du Congrês International de la Préparation des Minérais 
28, rue Arthur — Rozier 
Paris 19º, France 
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Vista do interior de uma das Fábricas "EFACEC" | | | | 
( Arroteia — S. Mamede de Infesta ) | | 
Distinguem-se nesta fotografia : 


— Um transformador de 15000 kV A/60 000 V (marca EFACEC) 
com regulação em carga, em curso de fabrico. 


— Um disjuntor de 220000 V (marca SOPREL), aguardando 


ensaio de rigidez. 


— Nota-se também um transformador de ensaios de 600 000 V | I ali 
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DESPOEIRAGEM 


FILTROS TECIDOS EM PEÇAS, PANOS, SACOS, MANGAS 


DE QUALQUER FIBRA NATURAL OU SINTÉTICA 


FANAF EL 


FÁBRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 


APART.: 9 
TELEF.: 983 OVAR 
TELEG.: FELTROS PORTUGAL 


Termómetros de mostrador para todos os fins industriais 
com a garantia NEGRETTI & ZAMBRA quanto a preci- 
são e robustez. 

À esquerda, um dos novos tipos “Standard” e, à direita, 
um dos modelos próprios para motores Diesel — apenas 
dois exemplos de uma larga gama abrangendo cente- 
nas de modelos. 

Grande variedade de modelos para entrega imediata 


na fábrica. 
NEGRETTI 


Sede: NEGRETTI & ZAMBRA LTD. 
Regent Street, London W.l 
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ESCRITÓRIO E FÁBRICA 


ESTRADA DE S. JOÃO 


& ZAMBRA 


Peça católogo em Português aos 

Representantes exclusivos em Portugal 

e Provincias Ultramarinas: 

MARTIN LESLIE & CIA. LDA, 

Campo das Cebolas, 43 3.º D.- LISBOA-2 
Tel. 32824 e 35234 


CG. D. U. 624.484 


TURBINAS A GÁS 


PELO ENG.º JOÃO CARLOS ADRIÃO DE SEQUEIRA 


Prof. do 1, 8. T. 


4.º lição proferida em Luanda e Lourenço Marques integrada 
no 1.º Curso de Férias da Universidade do Ultramar 


Princípio de funcionamento 


Um compressor C aspira o ar exterior comprime 
este ara 3 ou 5 atmosferas e envia o ar assim 
comprimido a uma câmara de Combustão B, ali 
é injectado óleo sob pressão que se queima com 
uma quantidade de ar em excesso determinado 
por forma a que a temperatura não exceda 600º 
centigrados nas instalações industriais, 800º a 
850º nos motores a jacto militares (Fig. 1). 


Fig. 1 


Os produtos da combustão a uma pressão 
elevada e a uma alta temperatura entram numa 
turbina T montada no mesmo veio da compres- 
são e do gerador eléctrico e fornecem na sua 
expansão um trabalho 30 “/y a 50º/o superior ao 
gasto na compressão e acabam por sair de novo 
para a atmosfera. É este excesso de potência que 
se vai encontrar na energia fornecida pelo ge- 


rador eléctrico montado no mesmo veio do com- 
pressor e das turbinas. 

Os gases dos escapes sairiam neste caso da 
turbina a 300º ou 400º sendo portanto suscepti- 
veis de aproveitamento principalmente por serem 
produtos da combustão efectuada com 60 “/o de 
excesso de ar e portanto ricos em oxigénio. 

O que se passaria neste conjunto seria repre- 
sentada no diagrama Pv e no diagrama entrópico 
pelos seguintes diagramas em que se admite que 
as compressões e expansões são adiabáticas. 


Fig. 2 


Temos assim de 1 a 6 compressão adiabática, 
de 6 a 9 combustão a pressão constante, de 9 a 
11 expansão adiabática na turbina e 11 a 1 escape 
para a atmosfera. 

Se em lugar da compressão ser adiabática ela 
fosse isotérmica a linha da compressão seria 
muito menos inclinada e o trabalho absorvido na 
compressão seria muito menor ficando portanto 
maior trabalho útil disponível como se pode ver 
no diagrama junto. 

Como a compressão isotérmica é impossível 
de se realizar adopta-se uma solução intermédia 
que é comprimir o ar em vários andares e arre- 
fece-lo nas passagens duns a outros andares da 
compressão. 


TEONIOCA 
589 


— — — Compressão adiabática 
— e Compressão isotérmica 


Fig. 3 


Por outro lado convém aproveitar o calor dos 
gases à sua saída da turbina e portanto usa-se 
um permutador do calor para aquecimento do ar 
necessário a combustão depois do último andar 
da compressão. 

O esquema das instalações passa a ser o das 
figuras juntas a que correspondem os diagramas 
entropicos e PV que se se seguem. 


Fig. 4 


Cj — Compressão da baixa pressão 
C, — Compressão da alta pressão 
Ki Ks K; — refrigerantes 

R —- Permutador recuperador 

T; — Turbina de baixa pressão 


T — » de alta pressão 
B, B; — Câmara de Combustão 
G — Gerador 


M — Motor de arranque 


Nas figuras temos a seguinte: 


Fig. 5 


(1-2) (3-4) Compressões dadas no compressor C, 
(2-3) (4-5) Arrefecimentos no arrefecedor K,; 

(5-6) Compressão no compressor C, 

(6-7) Arrefecimento no arrefecedor K; 

(7-8) (g-r0) (11-13) Compressão no compressor C; 
(8-9) (10-11) Arrefecimentos no arrefecedor K;» 


TRONTOA 
590 


(18-19) (13-14) Trocas de calor no Permutador R 
(14-15) Elevação de temperatura no queimador B; 
(15-16) Expansão na Turbina T, 

(16-17) Elevação de temperatura no queimador B, 
(17-18) Expansão na Turbina T, 

(19-1) Escape para a atmosfera 


cooper 


Fig. 6 


O emprego do permutador de calor R origina 
um aquecimento do ar necessário à combustão a 
300 e 350º sendo portanto só necessário empre- 
gar o combustível para levantar a temperatura 
a 600º. 

Os gases de escape vão baixar a sua tempera- 
tura e sair para a atmosfera a 150º e 200º a que 
se traduz numa franca melhoria de rendimento. 


DT =. 


TODOS 


Efe poa 
NR O 


a — compressor a baixa pressão 

b — refrigeração 

c — compressor a média pressão 

dd — compressor a alta pressão 

e — reaquecedor de ar 

f -- câmara de combustão à alta 
pressão 

g — turbina a gás a alta pressão 

h — câmara de combustão a baixa 
pressão 

ti — turbina a gás a baixa pressão 

k — motor de arranque 

! — redutor de velocidade 

m — alternador trifasado 

n — excitatriz 


cof 


É segundo um ciclo semelhante que trabalha a 
Central Suissa de Turbina a Gaz de Beznau 
que se vai passar a descrever. Como se vê na 
figura 6, temos 3 andares de compressão a c d in- 
tercaladas com 2 refrigerações intermédias (b) (b). 

O ar vai passar em seguida num premutador 
de calor e onde é aquecido com o calor dos 
gases do escape que saem atravez dele. 


F 


ig. 7 
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O ar vem em seguida entrar na câmara de 
combustão da alta pressão f e em seguida entra 


EROWM DOVEA ge TOati-d 
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na turbina da alta g que está montada no mesmo 
veio de compressor (d) do alternador trifásico 
(m) do motor de arranque (k) e da excitatriz (n). 

Os gases depois de trabalharem nestas turbi- 
nas vão passar na câmara de combustão de 
baixa pressão h de ali vão trabalhar na turbina de 
baixa i que tem montada no mesmo veio os 2 
compressores (a) e (c) e um outro motor de 
arranque k. 

Como se vê do gráfico junto (fig. 7 À) a recu- 
peração do calor permite atingir o rendimento 
de 30 “/o saindo os gases quentes a 180º. 

Chama-se a atenção para o númaro de rotações 
do grupo que chega a ser de perto de 3000 ro- 
tações por minuto. 

É igualmente para chamar desde já a atenção 
para a enorme bateria de refrigeradores que se 
foi obrigado a montar. 

A casa Brown Boveri constrói turbo-compres- 
sores denominados isotérmicos em que o ar sai 
a uma temperatura 35º a 40º superior à das 
águas de refrigeração. 

Nestes compressores não há refrigeração de- 
pois do 1.º andar por ser ainda a temperatura 
do ar muito baixa e não há refrigeração depois 
do úlmito andar por prejudicial, Esta refrigera- 


ção que de início era feita no interior do pró- 
prio compressor por circulação da água é agora 
sempre feita em refrigeradores colocados na vi- 
zinhança dos compressores para que o ganho 
que se fez na compressão não seja anulado pelas 
perdas nas condutas. 


dor 1 e se traçamos as curvas de rendimento 
do ciclo gás com compressão adiabática e com 
compressão isotérmica para temperaturas de gás 
de escape entre 550 a 700ºC mostra a figura 8 
em função de graus de compressão quais os 
rendimentos que é possível atingir nos dois 


Fig. 7c 


Em geral para pressão elevada usam-se 

9 andares de compressão com 7 andares de 
refrigeração 

6 ou 7 andares de compressão com 4 a 5 
andares de refrigeração. 


Se arbitrarmos como rendimento mecânico da 
compressão 0,8 da Turbina 0,85 e do alterna- 


— .— Compressão adiabádtica 


amipnacisanes ” isotérmica 
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5 10 15 20 
Relação de Compressão 
Fig. 8 


casos—e em que se vê que o rendimento pelo 
simples facto de se passar de compressão adia- 
bática para a isotérmica passa de 15 a 20 "fo. 
(Fig. 9) 

Como acabamos de constatar o ciclo da tur- 
bina a gás na sua forma mais simples não con- 
vém a produção económica da força motriz e 
só poder ser utilizada em casos especiais em que 
a economia não seja base absolutamente impres- 
cindível. 

Se porém examinarmos a temperatura a que 
os gases abandonam a turbina que são os da 
figura 9 vemos que se pode pensar claramente 
em melhorar o rendimento da turbina apro- 
veitando este calor no aquecimento do ar e 
reduzindo portanto o consumo de combustível. 
Este aquecimento tem de ser obtido em permu- 
tadores do calor e é óbvio que deve ser operado 
após as compressões e antes da câmara de Com- 
bustão com o esquema junto indica. 
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Rendimento 


Os gráficos juntos dão os melhorias dos ren- Tipos de compressores 
dimentos obtidos para uma saída dos gases das 


câmaras de combustão a 700º e para diferença À t Nas turbinas a gás usam-se 2 tipos de com- 


pressores; os compressores centrifugos munidos 
dos difusores com um único andar ou em mais 
andares em que a relação da compressão por 
andar varia entre 1,5 e 2,2 e que são do tipo 


700 
vulgar dos compressores centrifugos radiais como 
R o da figura junta em que o ar saído de um dos 
a andares é obrigado a dar uma volta de 180º 
y para entrar no andar seguinte como os das figuras 
a ia 12 e 13 e 14, e os compressores axiais hoje com 
300º emp. Gases 
entrada 
Ú 
ã 
É 
pd 100. 


E E) 8 1 
Rel. Comp. 


Fig. g 


da temperatura entre os gases à saída e o ar 
admitido entre 200º e 0º, gráficos que vêm mos- 
trar a vantagem enorme da recuperação do calor 
nos gases de escape. 


Fip. 10 


5 
Rei Comp. 
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Fig. 16 


largo emprego industrial e que devem o seu de- 
senvolvimento à sua aplicação nas turbinas a gás 
e que em princípio são uma turbina rodando em 
sentido contrario e transformando trabalho em 
pressão (Figs. 15 e 16). 


Os diagramas correspondentes a este tipo de 
compressores são em tudo identicos aos das tur- 
binas a vapor construindo-se estes compressores 
em 3 tipos: 


Compressores simétricos 
Compressores não simétricos 
Compressores Vortex 


O gráfico junto (Fig. 18) dá o número de anda- 
res que são necessários para obter vários graus 
de compressão nos vários tipos de compressores. 


30 


ba 
o 


Núm. Andares 


Rel. Comp. 
Fig. 18 


O quadro (fig. 19) estabelece a comparação en- 


Não simé 


, Simétrico | didi Vortex 
Número de rotações . 15.750 8750 5.000 
Número de andares .. 10 18 22 
Aumento de Pressão por 
andar . 1.149 | 1,08 1,062 
Rendimento cc. 82 87 | Bo 
Diâmetro Ext. . . 40 “/m 48 61 
E Diâmetro Base . .| 25ºm | 28 38 
E | Largura da Pá. .| 7,6%m | 10€/m| 12 c/m 
*: | Velocidade Perif.. «| 335m | 220m | 162m 
Velocidade Média. 160 m | 102m 61 m 
= Diâmetro Ext. . . «| 40 | 48 61 
T | Diâmetro Ext, « . «| 35 40 53 
5 Velocidade Perif.. 335 m/s| 220m/s| 162M/« 
É Largura da Pã. . . 2,5€/m| 3,7 €/m 4 “im 
A) 


144 m/s | 84 m/s 


Velocidade Média. 
Fig. 19 


tre os vários tipos de compressores para as mes- 
mas condições que são: 
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Relação de compressão 4k 
Pressão de admissão 1,033 
Temperatura de admissão 21º 


Pressão final 4k 
Temperatura final 126,º67 
Capacidade 679 mº/minuto 


Os diagramas dos vários tipos de compresso- 
res axiais são os indicados nos gráficos juntos 
(Fig. 17), em que o diagrama do andar simétrico 
está feito supondo 50 “/y de reacção. Fig. A,B,C. 


Diagrama da velocidade dum andar simétrico 


Fig. A 
C, — Velocidade entrada na directriz 
C. — Velocidade saída na directriz 
C, — Velocidade de entrada da roda motora 
C, — Velocidade de saída da roda motora 


Diagrama num andar assimétrico 
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Ao passo que na compressão centrifuga o ar 
tem que ser várias vezes desviado de 180º, no 
axial o ar tem um deslocamento muito mais sim- 
ples a partir dos 150.000 m” hora é de preferir 
sempre o compressor axial até esse valor, há que 
estudar caso por caso. O diagrama junto indica 
qual o débito mínimo a partir do qual se acon- 
selha a compressão axial. 


Examinados assim os vários tipos de compres- 
sores vamos passar às — Câmaras de combustão —, 

Umas vezes prevé-se uma câmara de combus- 
tão única como a da figura 21, em que o óleo 
combustivel é pulverizado por um queimador com 
ar em excesso para que a temperatura não exceda 
600º e em que a câmara de combustão é comple- 
tamente isolada para evitar perdas para o meio. 
O conjunto da turbina a gás torna o aspecto das 
figuras juntas em que se vê (fig. 22 e fig. 23). 


VIADUTO DE SACAVEM 
NA AUTA-ESTRADA DO NORTE 


Aspecto durante a execução dos ensaios de carga efectuados pelo Laboratório Nacional 


de Engenharia Civil 


dronstruçõe 
Y Técnicas, Ida 


SEDE: PRAÇA DO MUNICÍPIO, 13-3.º — LISBOA 
TELEFS. 366506 (4 LINHAS) 
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BOLETIM DAS VÁLVULAS CRANE DE CORREDIÇA EM BRONZE 


..apresentado a D46, com 
disco renovável em liga de 


Niquel e sede emaço inoxidável. 


Desenhada para trabalhar em condições 
realmente rigorosas e feita de acordo 
com as normas Britânicas 2060. O 
disco tipo cavilha, em liga de niquel 
e a sede em aço inoxidável tornam a 
vedação extremamente apertada ao mesmo 
tempo que permitem um controle rigo- 
roso do fluxo no estrangulamento. FEstes| 
dois metais diferentes proporcionam 
uma melhor resistência à corrosão e a 
matérias estranhas. Tanto o disco 
como a sede são renováveis. Pressões 
de trabalho (válvulas com as extremi- 
dades roscadas): até 250 libras (vapor 
500.ºF), 500 libras (água, óleo ou gás, 
frio, sem choque). 

Extremidades Hexagonais roscadas de 
acordo com as normas Britânicas-ou 
Americanas, Com flanges de acordo 
com as normas Britânicas, tabelas 'F! 
ou 'H' - ou com as normas Americanas. 
Dimensões: Roscadas 1/4'!' a 3!!-Com 
flanges l/2n a 5m 

Dimensões: 2!lm a 3n têm o castelo 
roscado. 


CRANE 


VÁLVULAS EM BRONZE, 
FERRO FUNDIDO E AÇO 


CRANE, LTD,, 15, 16 RED LION COURT, FLEET STREET, 
LONDON, E. C. 4 ENGLAND 


Agentes: R. E. B. E, L. 


Av. Joaquim António de Aguiar, 4 — LISBOA,1 
VA? 


Fig. 22 


O ar vindo do compressor de ar entra assim 
em todas estas câmaras de combustão dispostas 
em paralelo e segue depois para a turbina a gás. 

É este o tipo das câmaras de combustão empre- 
gados em todos os motores de aviação por ser 
mais difícil haver panes em todas as câmaras de 
combustão (Fig. 26). 


Turbina a gás — A queda entálpica que se obtém 
numa turbina a gás é sempre muito pequena. 
Por exemplo uma turbina alimentada com ar a 
4,812 aquecida a 600º tem à entrada no 1.º andar 
da alta uma entalpia de 213,3 k.cal por kilos de 
ar e à pressão do escape 1,803 é a 313º ainda 


it lno; + [oad 


Fig. 23 


1 — Gerador eléctrico 
2 — Turbina alta pressão 
3 — Câmara de combustão 


Outras vezes em lugar duma só câmara de 
combustão empregam-se muitas câmaras em para- 
lelo, algumas vezes 12 ou mais conforme se vê 
na gravura junta (Fig. 24) tendo cada câmara de 
combustão com um diâmetro de 112 mm e em 
comprimento de 900 mm sendo o envólucro ex- 
terior em aço de carbone com um dispositivo 
em fole para a dilatação. A camisa interior dentro 
da qual se dá a combustão é em aço cromo 
níquel soldado numa das extremidades e livre na 
outra para dilatação. 

O óleo é injectado no topo da câmara de com- 
bustão e o ar primário é ali igualmente injectado, 
o ar secundário entra por uma série de orifícios 
no tubo interior e é portanto aquecido ao passar 
no intervalo dentro do tubo (Fig. 25). 


4 — Turbina baixa pressão 
5 — Compressor de ar 
6 — Motor de arranque 


7 — Permutador de calor 
8 — Conducta de aspiração 
9 — Conducta de escape. 


tem 144 calorias o que quer dizer que a queda 
entálpica em toda a turbina que trabalhasse nestas 
circunstâncias seria 


213,3 — 144 = 69,3 calorias 
Como em geral estas turbinas são divididas 


TEONICA 
597 


Fig. 25 


Fig 4 Dispensação doe pump 
sentes do anoror de tres bemjaso 
[5 Co ecremdrios fome & qu 
eguveda, ec Isso da qcutemumi- 
Sarto dr ademiador. com mw cbimia- 
vo de comido mo eópigu po 
river de suriuas db direma 


Fig. 27 


Diagrama de um andar 


em alta e baixa por causa das temperaturas rei- 
nantes e como o número de rodas motoras oscila 
entre 10 e 11 na totalidade vem uma queda entál- 
pica por andar da ordem das 4,2 calorias kilo. 
Quase todas as turbinas de gás são de reacção 
com coeficiente de reacção variável entre 0,45 
e 0,65 crescentes para a saída dos gases. 
Teremos assim para um andar utilizando uma 
queda de 4,2 calorias uma queda de 
Ah = 0,45 x 4,2 = 1,89 na directriz 
Ah” = 0,55 X 4,2 = 2,31 na roda motor 


a velocidade de escoamento dos gases na 1.º roda 
directora vem portanto a ser 


Ca = 91,53 2,31 = 140 m 
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a velocidade real da saída seria devido aos atritos 
na directora 


C,=0,97 < 140 = 135m. 


O ângulo de saída e adoptada tanto nas direc- 
trises como na roda motora é de 20º. O número 
de rotações de turbina é em geral de 3000 que 
com um diâmetro médio de 0,925 m dá uma 
velocidade 

u == 145 m/segundo 


valor que me permite fezer o diagrama de entrada. 

Na roda motora a gás vai aumentar a veloci- 
dades ao longo da pá devida a queda de 1,89 
calorias 


Wa = 0,97 x 91,53 V 1,89 + 


E + di 
426 x 2x 9,8 


W: = 0,97 < 91,53 V/1,89 + 03 = 134 m 


Valores que nos mostram ser as pás muito 
direitas tanto da directriz como das rodas motorra, 
visto que a directriz varia os ângulos entre [2 e a 
e a roda motora entre [4 e as. 

Como o volume de gases é muito grande a 
altura das pás é grande as diferenças de veloci- 
dade v são grandes o que faz com que as pás 
sejam todas torsas como se vê na figura 28. 

O número de andares é sempre pequeno nestas 
turbinas. 


Fig. 28 


a — Compressão axial 

b — Câmara de combustão 
c — Queimador 

d — Entrada do ar ciclônico 
e — Envólucro arrefecedor 
f — Palhetas da turbina 

g — Turbina a gás 

h — Válvula de segurança 
t — Alternador 

k — Motor de arranque 


Fig. 28a 


Materiais para as pás das rodas motrizes 


Apesar de propositadamente se ter feito a com- 
bustão com um grande excesso do ar para que 
a temperatura na admissão não excedesse 600º 
ou 800º, como o número de rotações é sempre 
elevada 3000 r. p.m. vão as pás das turbinas estar 
ao mesmo tempo submetida a altas temperaturas 
e a grandes esforços de tracção devida à força 
centrifuga tendo portanto de se ter em conta a 
fluage tanto mais de temer quanto maior for 
o número de rotações como no caso presente. 

A resistência à fluagem deve ser tal que nunca 
haja o perigo do material das pás alcançar a 
fluagem terciária durante a vida das máquinas 
porque isso representaria a desagregação com- 
pleta da roda motriz. 

Em virtude deste facto os aços correntes empre- 
gados na confecção das pás das turbinas a vapor 
não servem. Além da resistência à fluagem o 
material das pás deve apresentar as seguintes 
caracteristicas : 


— Estabilidade da estructura a temperaturas 
elevadas. 

— Boas resistências às vibrações e portanto a 
esforços alternados. 

— Facilidades de fabrico. 


Os materiais ou são aços ligas ou ligas metá- 
licas em que o ferro entra em quantidades mí- 
nimas. 

Para os aços usam-se aditivos que são: 


Cobalto — 12 a 14 º/o que dá uma grande resis- 
tência à oxidação em escamas. 

Diminui a corrosão a altas temperaturas. 

Silício — 3 */) aumenta a resistência a corrosão 
em escamas «scaling». 


Os aços inoxidáveis cromo-níquel como os 
empregados pela Brown Boveri. 


C —- 0,2 
Si — 1,2 
Nec — 8 

Cr — 18 
W —1 

Mg — 0,6 
fe — 71 


são de empregar até 600º q, da ordem dos 750 
a 800 kem”, 
Os ingleses usam para esse efeito ligas for- 
jáveis de 
80 “/a — Níquel 
20 "/y — Cromio. 


Um outro material muito usado é o «Nimonic 
80» que resiste bem a 750º de temperatura com 
a seguinte composição 


C — 0,04 
Mn — 0,56 
Si — 0,48 
Fe — 0,44 
Cr — 21,18 
Ti — 2,44 
Al — 0,63 
Ni — 7423 


Este material tem uma grande resistência à 
fluagem quando tratada da seguinte maneira: 


9 horas a 1080º C 
arrefecimento ao ar reaquecimento 
16 horas a 700º €. 
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Convém empregar este material por forma que 
o alongamento por fluagem não ultrapasse durante 
a vida das pás da turbina 


0,1 0 


Os ensaios com este material indicaram que 
este alongamento é obtido em conformidade com 
o gráfico da fig. 29. 


1009 5000 10000 50000 100000 horas 


Fig. 29 


Os dados sobre experiências prolongadas com 
materiais resistentes a altas temperaturas mos- 
tram que uma vida de 10000 horas é compatível 
com uma carga 945 kcm? para o «Nimanic 80» 
até 700º. 

Em motores a jacto de avião pode-se ir para 
maiores temperaturas e cargas mais elevadas 
diminuindo as horas de vida atribuíveis às pás. 


Sobe-se assim com a fadiga a valores até 1800 
kcm? nos motores a jacto diminuindo as horas 
possíveis de trabalho. 

A maioria das pás das rodas motoras são 
obtidas das peças forjadas em bruto e depois 
acabadas cuidadosamente. 

Um outro grupo de pás é feita em ligas me- 
tálicas duras que dificilmente se deixam traba- 
lhar, recorre-se então à moldação em cera per- 
dida que permite 


— Limites de precisão de 0,05 mm 

— Custo da produção reduzido 

— Tendo como limitação os pesos do Metal 
que em geral não pode exceder 6,8 kg. 


Motor a jacto 


É uma turbina a gás em que a potência desen- 
volvida pela turbina é absoluamente igual à 
necessária ao compressor do ar. 

Vão assim os gases quentes abandonar a tur- 
bina através da duse de saída com uma grande 
velocidade desenvolvendo uma potência 


F Vi 


Sendo F a força e V; a velocidade de saída. Se 


MATERIAIS PARA PALHETAS 


| Mn Si Ni Co Cr | Mo | W | Cb | Ti Al | Forro 

E— 4% E. E mei 0,78, 0,341 41,3 |21,7 117,7 e a E 2,55 0,55! 15,3 
Refractaloy 26 | 0,04 0,8 [0,93 37 | 20 | 18 | 3,2 | —- | —s | 2,52 0,24] 17 

» 70 .|011/1,97/0,25/201 |30,1 |20,2 | 83 |38 | — | — | — | 152 
Discoloy 26 .. | a 0,7 |0,4 |25,2 sima 8,3 sao | — a 233, 0,42 | Restante 
S 816 0,4 | 0,58  0,54|20,23|47,7 |19,5 3,93/3,45|406| — | — | 2,95 
Niamonic 80 0,04 0,56. 0,48/74,23] — |21,18] — | — | — |2,44/0,63| 0,38 
Incanel W . 0,031 0,53] 0,52 /7482] = |T4A | = | = | —-— | 277] 068] 617 
Rolled Vitalium .|0,22/ 0,73 0,57| 2,92 61,4 |27,56| 59 | — cms — em 0,50 
172 W.....|0,46|0,58| 0,43|19,16| — |13,19 0,52. 248) — | — | — | Restante 
Tinken 16-25-6 | 0,08/1,34| 0,31/23,15| — |17,95| 6,79] — Ea E 51 
Hasteloy B . 0,05 0,59 0,19 | 65,10| — o | 28,683] su | a | quai) 47 
Rex 78 «su 0,08 sa 18 — | — — | 0,6 — | Restante 
Tinidur T. ../0,15/0,6 10,7 | 30 | — a 15 — | — [1,7 | — |(Ca 4.0) 
NISG cs... Dad | LAS an) ançaa casa [a 66| 3,02 | 2,10 108) — | — A 

Fig. 30 
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para se produzir este jacto se gastou uma ener- 


gia P o rendimento da propulsão 
a = FM 
P 


Se o jacto sair do aeroplano com uma velo- 
cidade v; e se o aeroplano se deslocar com uma 
velocidade v e se através do jacto se escoarem 
Wa quilos será 


e a energia perdida no jacto será 


P, = fiz (vi — v)* 
2 8 
e o rendimento 
O ui FEV 
P, + EV 


por ser P;4+FV a energia total. Substituindo 
valores vem 


Ma (vi— v) v 
g | 


A ES EEE 


e (vi — v)* + Ea (vj— v) v 
2g 8 


OU 
Bl) 


(vy—v)? + 2(vi—v)v 


2v 2v 
à DESSE A E ES 
vVi— v) + 2v v+ovi 
Da expressão 
W 
E = = (y; —v) 
8 
tira-se 
E: 
ViI— V= a 
Wa 
4 D sus: - Fy 
v Wav 
Vi F 
Do mid do 
V Wa V 


Expressões que permitem determinar a Força 
do jacto e rendimento da propulsão para várias 
velocidades. 


Motores de propulsão por jacto 


Os dados e gráficos que a seguir vamos pro- 
jectar referem-se a motores Americanos e Ingleses 
de tipo não secreto. Há em todos os motores a 
jacto modernos duas tendências claras: 


— Emprego de compressor axial em lugar de 
compressor centrifugo 


— Aumentar a taxa de compressão. 


— À custa da vida do motor aumentar a tem- 
peratura de combustão a 815º C na admissão das 
turbinas valor este que junta ao da taxa de com- 
pressão originam um decréscimo no consumo em 
kilos combustível por hora e kilo de propulsão, consumo 
ao nível do mar que tem tendência a descer 
a 0,7 k (Fig: 31). 


g propulsão 


kg combustão/hora/k 


O peso de motor por kilo de propulsão (Fig. 
32) tem uma tendência a diminuir tendo passado 
0,55 kilos para 0,3 k do peso por kilo propulsão. 

ass 


050 


040 


030 
1943 1947 1949 1951 1953 1955 
Fig. 32 
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Um outro gráfico interessante é o coeficiente 
de propulsão por unidade de superfície frontal; 
este gráfico mostra a tendência clara do au- 
mento e este muito mais acentuado nos motores 
de compressão axial de que nos motores de com- 
pressão centrífuga passadas nas primeiras de 
0,196 k a 0,296 k e nas segundas de 0,098 k 
a 0,196 k (Fig. 33). 


0147 


0.098 


1943 1947 1949 1951 1953 1955 


Fig. 33 


Como já vimos a impulsão é igual à massa, 
multiplicada pela velocidade de escape. Para 
aumentar a força da propulsão é necessário ou 
aumentar a massa ou aumentar a velocidade de 
escapes. 

Um dos processos usados é injecção de água 
para refrigerar os gases antes da entrada na 
turbina, a água permite queimar mais combustível 
sem exceder a temperatura tolerável aumentando 
a massa dos gases. Pode por esta maneira aumen- 
tar-se a força na propulsão em 20 */', 

Um outro processo que se usa é a combustão 
adicional no tubo de escade após a turbina a que 
pode originar aumento da potência da ordem dos 
40 */o sendo muito usado este processo em aviões 
de caça, 

As primeiras provas levados a efeito com 


motor de jacto mostraram que à proporção que * 


a altitude aumentava com diminuição da tempe- 
ratura e da densidade do ar, mais difícil se tor- 
nava manter uma boa combustão diminuindo o 
rendimento da câmara de combustão até extin- 
guir-se as chamas a 9000 m de altitude. 

Dum modo geral a estabilidade da chama 
melhora quando decresce a velocidade dos gases 
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na câmara de combustão o que leva a aumentar 
o tamanho das câmaras de combustão e portanto 
o peso para motores a jacto destinados a grandes 
altitudes (Fig. 34). 


Jacto 


Turbina 
Fig. 34 


A compressão é parcialmente obtida pela 
máxima velocidade do avião que dando origem 
a uma entrada de ar origina uma pressão inicial 
que aumenta com a velocidade do avião. O dia- 
grama entrópico do motor a jacto é assim dado 
pelo seguinte: 

A energia que se obtém devido a compressão 
inicial do ar é encontrada no jacto (Fig. 35). 


E ner.Tur bina 
CompsHélice E CompsHélice 


Jacto 


Fig. 36 


Outras vezes usa-se o sistema misto de jacto 
e propulsão por hélice como vai indicado na 
figuras 36 e 37. Neste caso a turbina a gás tem 
de ter uma potência igualàã energia gasta no com- 
pressor mais a gasta na hélice sendo portanto 
a energia do jacto mais pequena. Igualmente em 
voo a entrada do ar aumenta a força do jacto. 

Já tivemos ocasião de chamar a atenção para 


velocidades -=»m————————— Em Vôo 


Fig 37 


o facto dos produtos de combustão das turbinas 
de gás sairem a uma temperatura elevada e 
serem ricos em oxigênio em virtude do excesso 
de ar a que se foi obrigado para reduzir a tem- 
peratura dos produtos de combustão a 600º. 

Isto torna possível empregar os gases de 
escape das turbinas a gás como ar de combus- 
tão em caldeiras conduzindo a um aproveita- 
mento total do calor nelas contido. Desta apli- 
cação resulta que o rendimento da produção da 
energia será 


N 


u=860 —————— 
O ot — Qr 


em que 


N -— A potência útil da Turbina a gás kW 

Qtot— Calor total fornecido K.cal/hora 

Or — Calor que seria fornecido a uma for- 
nalha independente da mesma potên- 
cia K.cal/hora 

860 — Factor de conversão KWh/K. cal hora. 


As figuras juntas dão uma ideia da forma 
mais simples de uma tal combinação. 

É claro que nestas aplicações não há qualquer 
interesse em fazer as compressões com arrefe- 
cimentos finais. 


O rendimento duma instalação deste gênero é 
sempre melhor do que os duma Central Térmica 
independente e varia entre 40 e 80º/0. 

As considerações que acabamos de fazer em 
nada são alteradas se o ar comprimido em lugar 
de ser aquecido numa câmara de combustão o 
for por intermédio de vapor e se a turbina pas- 
sar a ser de ar quente em lugar de gás. 


Perm. 


Fig. 39 


O aproveitamento dos gases de escape da tur- 
bina a gás para produzir vapor numa caldeira de 
recuperação traz consigo a dificuldade de não 
poder arrefecer os gases a uma temperatura sufi.. 
cientemente baixa por terem estes ao sair a uma 
temperatura sempre bastante superior à da água 
da alimentação. 


Fig. 40 


TV; — Turbinas a vapor 


TVa a » » x» 

C | — Compressor 

B  — Câmara combustão 

TG — Turbina de gás 

R  — Permutador 

P | -— Bomba alimentação 

G | — Geradores eléctricos 


TÉONIOA 
603 


O débito da água determinada pelas possibili- 
dades da evaporação é neste caso insuficiente 
para permitir um arrefecimento dos gases no 
aquecedor de água. O problema agrava a sua 
complicação à medida que a pressão do vapor 
vai subindo e que a temperatura dos gases à 
entrada na caldeira vai baixando o que contra 
indica o emprego da caldeira de recuperação num 
cíclo a gás. 

Há portanto necessidade de encarar o problema 
de outra maneira, 

Determina-se o débito de água necessária para 
que os gases de escape saiam à temperatura con- 
veniente a um bom rendimento. A saída do aque- 
cedor da água R vamos obter ou água quente 
ou mesmo uma mistura de água a vapor. O Per- 
mutador R não é uma caldeira de recuperação 
mas sim um aquecedor da água da alimentação 
em que este aquecimento pode ser bastante ele- 
vado pois se pode nele obter uma mistura de 
agua a vapor. 

No caso de se dar só aquecimento da água 
supondo que a água vinda dos condensadores 
chega a 25º e que se pretende arrefecer os gases 
a 40º vem que a quantidade da água que se pode 
aquecer por kilo de gás é 


g (de — 25) = (t — 40) Cp 


sem que +e é o calor do líquido e t a temperatura 
da saída dos gases e Cp o calor específico dos 
mesmos. 

No caso de se querer produzir já vapor satu- 
rado no permutador o que não seria de aconse- 
lhar teriamos 


g (Je — 25) = (t' — 40) Cp 
g(A—de)=(t—t)€, 


em que » é o calor total de vaporização e tí a 
temperatura dos gases no ponto em que começava 
a vaporização. 

A câmara de combustão da turbina a gás é 
agora a caldeira onde se ultima a vaporização 
e onde se vai dar o sobreaquecimento com a 
vantagem de que os produtos da combustão vão 
servir na turbina a gaz e sair finalmente a 40º 
evitando-se as perdas em calor sensível pela cha- 
miné. Como se acaba de ver a turbina a gás apa- 
rece aqui com um complemento da instalação de 
vapor melhorando-se o rendimento do ciclo de 
vapor porque da combinação caldeira-câmara da 
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combustão resulta uma caldeira de alto rendi- 
mento. As curvas juntas dão os rendimentos das 
ciclaã mistas gás e vapor que como se vê atin- 
gem 0,39. 

É fácil numa pequena aplicação examinar os 
resultados do sistema. 


Pressão da caldeira 40 kcmº 
Je a esta pressão 301 cal kilo 
À » » » 662 cal kilo 


Supondo o compressor com uma relação de 
compressão de 5, e uma compressão adiabática 
encontra-se pela fórmula 


g (he — 25) = (t — 40) Cp 
que 


g = 0,327 k vapor/kilo de gás 


Rendimento 


o gráfico junto indica o peso do ar em kilos por 
kW hora produzido e por ele se tira que nas con- 
dições estabelecidas a quantidade de gás seria no 
caso indicado de 39,4 k por kW hora e que por- 
tanto a quantidade de água que se poderia aque- 
cer seria de 


39,4 x 0,327 = 12,88 k de água por kW hora gás. 


Como o consumo de uma turbina a vapor a 
65k e 500º é em geral de 3,36 k/kW hora a pro- 
dução de vapor seria de 


Elia 


— 3,835 kWh vapor/k W h gás. 
irá por/ gás 


A quantidade de calor necessária a esta pro- 
dução de vapor com água já aquecida a 300º é de 


12,88 (817,2 — 301) = 6.660 cal 
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FABRICAÇÃO. SEMPRE RECENTE 
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| FLUIDOTERMUS 


Aquecedores-acumuladores industriais 


Dad 
a 


A disposição da Indústria maclonal 
Processo patenteado francês de aquecimento por fluido especial transmissor de energia térmica 


—- Temperatura até 300º C 

— Não têm pressão (Não sujeitos a prova) 

— Não têm incrustações nem oxidações 

— Não têm purgadores 

— Não têm panes nem periodos de paregem. 


Com um pequeno permutador a mesma instalação pode produzir vapor 


6 unidades para uma produção de 6,000.000 k cal h. foram já encomendadas pela indústria petrolífera para Portugal. 
Fabricação inteiramente nacional 


VASCO PESSOA, LDA. 


Rua da Boa Vista, 63 —- LISBOA 


IMPERMEABILIZAÇÕES 


TORNA TUDO ESTANQUE n 
Impermeabilização e vedação perfeita de: 


Terraços, Empenas, Algerozes, Clarabóias, 
e BIO Paredes húmidas, Coberturas, Fendas, 
Juntas de dilatação, etc., etc. 


Fácil aplicação = Elasticidade excepcional 
Não endurece nem apodrece 


FLUIDO E MASTIC 


RODFATEX 


Reforça a acção impermeabilizante do «Rubio», dando às construções um aspecto agra- 
dável — Protege todos os materiais contra a humidade, a ferrugem e a corrosão 


À base de borracha + alumínio 


RUBIO E ROOFATEX dois produtos que se completam 


Representante Exclusivo: (olundo TFemandes 


Rua Palmira, 46-A-B-C — LISBOA-1 — Telefs. 83 0163-8476 39 — End. teleg. «ORFER» 
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